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ABSTRACT 

We self-consistently derive the magnetic energy and relative magnetic helic- 
ity budgets of a three-dimensional linear force-free magnetic structure rooted in 
a lower boundary plane. For the potential magnetic energy we derive a gen- 
eral expression that gives results practically equivalent to those of the magnetic 
Virial theorem. All magnetic energy and helicity budgets are formulated in terms 
of surface integrals applied to the lower boundary, thus avoiding computation- 
ally intensive three-dimensional magnetic field extrapolations. We analytically 
and numerically connect our derivations with classical expressions for the mag- 
netic energy and helicity, thus presenting a so-far lacking unified treatment of 
the energy /helicity budgets in the constant-alpha approximation. Applying our 
derivations to photospheric vector magnetograms of an eruptive and a nonerup- 
tive solar active regions, we find that the most profound quantitative difference 
between these regions lies in the estimated free magnetic energy and relative 
magnetic helicity budgets. If this result is verified with a large number of active 
regions, it will advance our understanding of solar eruptive phenomena. We also 
find that the constant-alpha approximation gives rise to large uncertainties in 
the calculation of the free magnetic energy and the relative magnetic helicity. 
Therefore, care must be exercised when this approximation is applied to photo- 
spheric magnetic field observations. Despite its shortcomings, the constant-alpha 
approximation is adopted here because this study will form the basis of a com- 
prehensive nonlinear force-free description of the energetics and helicity in the 
active-region solar corona, which is our ultimate objective. 

Subject headings: MHD - - Sun: atmosphere - - Sun: chromosphere - - Sun: 
corona — Sun: magnetic fields — Sun: photosphere 
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1. Introduction 

The magnetic origin of solar eruptions has been established over the past several decades 
of solar research. Most eruptions originate in active regions that are, in general, closed 
magnetic structures rooted in the solar photosphere. Magnetized plasma motions in the 
solar atmosphere prevent a magnetic structure from attaining a minimum-energy, current- 
free state. Excess magnetic energy in active regions is manifested by the appearance of 
electric currents (Leka et al. 1996). Eruptive and noneruptive manifestations, such as 
coronal mass ejections (CMEs) and confined solar flares, respectively, must be fueled from 
this reservoir of free magnetic energy that is thought to be released in intermittent episodes 
of magnetic reconnection. 

A popular, early, view of a magnetic energy release event involved a nonpotential pre- 
event state relaxing into a potential, or nearly potential, post-event state of the magnetic 
configuration. The excess (nonpotential) magnetic energy was thought to be released during 
the relaxation. Seminal works on magnetic helicity (Woltjer 1958; Calugareanu 1961; Berger 
& Field 1984; Finn & Antonsen 1985; Berger 1985; 1988; Moffatt & Ricca 1992 and others), 
however, demonstrated that this view is incomplete or even misleading: magnetic helicity 
relates to the linkage of a magnetic structure (twist, torsion, and writhe), which is globally 
invariant even under resistive processes, such as magnetic reconnection. Helicity is present 
wherever electric currents are present. Therefore, an isolated helical magnetic structure 
cannot relax to a potential state unless its magnetic helicity is bodily removed from it. 
This provides a plausible interpretation for CMEs (Low 1994; Rust 1994a; 1994b), provided 
that the magnetic helicity in the erupting structures is not transferred to other parts of 
the solar atmosphere along preexisting or reconnected magnetic field lines. Indeed, eruptive 
activity seems to be necessary for the Sun where differential rotation and subsurface dynamo 
continuously generate helicity in the two solar hemispheres (Berger & Ruzmaikin 2000), with 
a statistical hemispheric segregation exhibited by magnetic structures of opposite senses of 
helicity (Pevtsov, Canfield, & Metcalf 1995). In confined events the magnetic configuration 
cannot relax to the potential state but, at best, it may relax to the lowest possible energy state 
that preserves the pre-event amount of magnetic helicity. This is known to be a constant- 
alpha, linear force- free (LFF) state (Woltjer 1958; Taylor 1974; 1986). Several works rely on 
the Woltjer- Taylor theorem, although controversy remains over its applicability to the Sun 
(Kusano et al. 1994, but also Antiochos & DeVore 1999). 

Whether magnetic helicity per ce is important for solar eruptions is also a subject of 
debate (e.g. Rust 2003; Rust & LaBonte 2005, but also Phillips, MacNeice, & Antiochos 
2005). Regardless, however, knowledge of the magnetic helicity is essential for a complete 
assessment of the magnetic complexity present in the solar atmosphere. Berger & Field 
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(1984) and Finn & Antonsen (1985) derived a gauge-invariant definition of magnetic helicity 
applying to open and multiply connected volumes such as the ones assumed for the solar 
atmosphere. The resulting relative magnetic helicity subtracts the helicity of the reference 
(potential) field so a nonzero value implies by definition the presence of free magnetic energy 
in the configuration. The relative magnetic helicity H m has two equivalent forms in the 
above works, namely 



where B p and A p are the potential magnetic field and its generating vector potential, re- 
spectively, and B, A are the respective quantities of the nonpotential field. The integration 
refers to the open volume V that contains the part of the magnetic structure extending above 
a lower boundary. To derive equation (JTJ, both Berger & Field (1984) and Finn & Antonsen 
(1985) assumed nonlinear force-free (NLFF) magnetic fields. The force-free approximation is 
probably necessary for helicity calculations because only in this case one obtains some knowl- 
edge of the magnetic field vector and the generating vector potential required to evaluate 
equation ([!]). 

Equation ([!]) cannot be evaluated in the active-region atmosphere, however, because 
the magnetic field vector is unknown above the lower boundary, be it the photosphere or 
the low chromosphere. Currently, active-region magnetic fields can only be measured in this 
boundary, so the only way to evaluate the relative magnetic helicity through equation (JT]) 
is by force-free (preferably NLFF) field extrapolation into the active-region corona using 
the measured magnetic fields as the required boundary condition. However, the NLFF 
extrapolation of observed solar magnetic fields remains an active research area where even the 
most successful of the existing techniques (Schrijver et al. 2006 and references therein) are too 
slow to fully exploit the spatial resolution of modern (let alone, future) magnetographs. The 
NLFF approximation should always be pursued given that the LFF approximation is almost 
certainly an oversimplification for most active-region fields and it can even be misleading 
in several cases. Even the NLFF approximation is most likely invalid in the photosphere 
(Georgoulis & LaBonte 2004), although it may hold in and above the chromosphere (Metcalf 
et al. [1995]; see, however, Socas-Navarro [2005]). 

Even in case equation (jTJ) is evaluated, however, it does not establish a link between 
the relative magnetic helicity and the magnetic free energy of the studied configuration. 
In addition, it might be risky to evaluate a volume integral of extrapolated fields at large 
heights above the boundary because numerical effects might settle in and affect the result. 
One, therefore, envisions a convenient surface-integral representation of the relative magnetic 
helicity that might alleviate the need for full-fledged three-dimensional extrapolations. To 
our knowledge, this has been attempted only in the LFF approximation following either the 
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theoretical analysis of Berger (1985) or the "twist" helicity of Moffatt & Ricca (1992). In the 
first case (Demoulin et al. 2002; Green et al. 2002), the employed formula for the relative 
magnetic helicity is 

where a is the unique, representative value of the force- free parameter and h UuVm is the 
Fourier amplitude of the measured normal magnetic field for the harmonic (ui,v m ) in a 
two-dimensional Fourier space with linear dimensions n x , n y . Equation (j2J) is a linearized 
version of the actual formula of Berger (1985). Linearization helps avoid \H m \ — > oo when 
\a\ — > (2tc/L), where L is the linear size of the computational domain. Detailed discussions 
and an extension of Berger's (1985) analysis will be given in §4.2 and Appendix B. In the 
second case (Regnier, Amari, & Canfield 2005) the magnetic helicity is approximated by the 
"twist" helicity of a semi-circular, constant-alpha magnetic flux tube, i.e., 

H m = \aL& , (3) 

o 

where L is the characteristic footpoint separation length of the tube and $ is the magnetic 
flux carried by the tube. Though useful, equations (T5]) and ([3]) also lack a much wanted link 
between the relative magnetic helicity and the free magnetic energy in the LFF magnetic 
structure that would enable a complete, self-consistent, description of its energetics. More- 
over, it is not clear how to generalize equations §Z§ and ([3]) into a NLFF calculation that, as 
should be always kept in mind, must be the ultimate objective of the calculation. 

The above difficulties and lack of information in the calculation of the total relative 
magnetic helicity prompted alternative lines of research. The lower boundary of a closed 
magnetic structure, where all magnetic field lines are supposed to be rooted, acts as the driver 
of the evolution in the structure either via boundary flows or via the injection of additional 
structure through it. Therefore, magnetic helicity can either be transported to and from the 
structure through this boundary or it can be generated by flows on the boundarjQ. Based on 
these principles, Berger & Field (1984) derived a surface- integral expression for the temporal 
variation (dH m /dt) of the relative magnetic helicity in a magnetic configuration. Besides its 
dependence on magnetic field vectors and vector potentials, as in equation ([T]), (dH m /dt) 
depends on the boundary flows. The advantage of the Berger & Field (1984) expression for 
(dH m /dt) is that it does not explicitly require force-free fields. The calculation of (dH m /dt) 
has been attempted by numerous authors over the past few years (see, e.g., Nindos 2006; 



■"■In the Sun, generation of helicity above the photosphere automatically implies the generation of an equal 
and opposite amount of helicity below the photosphere, to ensure a zero net helicity. 
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LaBonte, Georgoulis, & Rust 2007, and references therein), although it suffers from the lack 
of a reference value, namely the total relative magnetic helicity. The total relative helicity is 
a focus of this work, so the formula of Berger & Field (1984) will not be discussed further. 

This study is the first of a series of studies that perform a self-consistent calculation of 
the total magnetic energy and relative magnetic helicity in a closed magnetic configuration. 
We devise a practical way to calculate magnetic energies and helicities from solar active- 
region vector magnetograms provided that the observed magnetic configuration is isolated 
and flux-balanced on the presumed "plane" of the observations. The final expressions for 
the magnetic energy and helicity are derived in the form of surface, rather than volume, 
integrals. We always assume that the studied magnetic field configuration is in a force- 
free equilibrium. In this study, we provide the analytical foundation of a constant-alpha, 
LFF, energy-helicity calculation. A NLFF generalization of the energy-helicity equations 
will be the subject of a later study. The LFF analysis of this work relies on the energy- 
helicity formula of Berger (1988) evaluated via an application and extension of Berger's 
(1985) analysis for the magnetic energies and the relative magnetic helicity. Our objective in 
this work is the derivation of practical LFF energy and helicity equations that can be readily 
applied to solar vector magnetogram data. The magnetic energy budgets for a constant-alpha 
magnetic configuration are discussed in §2. The LFF energy-helicity formula is discussed in 
§3. The relative magnetic helicity is derived both as a volume and as a surface integral in 
§4. In §5 we apply our LFF analysis to vector magnetograms of two solar active regions and 
in §6 we summarize and discuss our analysis and results. 



2. Gauge-invariant definitions and the magnetic energy equation 

2.1. Magnetic field and the vector potential 

Assuming planar geometry, consider a magnetic structure B extending in the half space 
z > above a lower boundary S (z — 0). Let an open volume V of the half space z > 
include the structure and extend to infinity with its only boundary being the surface S. If 
S is not a flux (magnetic) surface, i.e. if B • z\$ ^ 0, then the configuration is analogous 
to a solar magnetic structure rooted in a small (assumed planar) part S of the photosphere 
and extending to infinity above S. Here z is the unit vector along the z-axis of a Cartesian 
coordinate system with an arbitrary origin on S. In the absence of plasma, this magnetic 
configuration can only be the vacuum, current-free magnetic field B p if the configuration is 
isolated (not interacting with other configurations) and flux-balanced on S. The presence of 
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plasma dictates a current- carrying magnetic structure B c , such that 

B = B P + B C . (4) 

The divergence-free properties of B and B p together with equation (j3j) ensure that B c is 
also divergence-free. As a result, we can define generating vector potentials A p , A, and A c 
for Bp, B, and B c , respectively. By definition, 

V x A p = B p , (5a) 

V x A = B . (5b) 
In addition, the Coulomb gauge is adopted for both A p and A to provide 

V ■ A p = , (6a) 

V • A = . (6b) 

Since S is not a flux surface, however, both the definition of B p and some topological 
properties of the field, most notably the ones present in its magnetic helicity integral, are 
not unique (gauge-invariant) and hence lack a physical meaning (e.g. Dixon et al. 1989; 
Berger 1999). Berger (1988) addressed the problem by providing gauge conditions for A and 
A p such that both B p and the magnetic helicity can be uniquely defined. These conditions 
are A p • n\g V = and A p x h\g V = A x n| 9v and were formulated for a volume V bounded 
by a surface dV, where n is the unit vector normal to dV and oriented outward from V. 
If V extends to infinity, Berger (1988) stresses that A p and A must additionally vanish at 
infinity. This restricts the above gauge conditions to the lower boundary S, so that 

A p -z| 5 = , (7a) 

A p x z\ s = A x z\ s , (7b) 

where z = — n. From equation (j4]) and the conditions of equations (5)- (7) we can derive 
additional conditions for the vector potential A c . Writing equation (j3j) in terms of vector 
potentials, we obtain A = A p + A c + V0, where <fi is an arbitrary scalar. Choosing the 
gauge such that V0 = 0, one obtains 

A = A P + A C , (8) 
If one now takes the dot and cross products of equation (jSJ) with z, then one obtains 

A c • z\ s = A • z| s , (9a) 
A c x z| 5 = , (9b) 
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where we have used equations (7a) and (7b) to reach equations (9a) and (9b), respectively. 

For a volume V bounded by a surface dV, another condition for the uniqueness (gauge- 
invariance) of B p and the magnetic helicity is that B and B p share the same normal compo- 
nent on dV, i.e. B • n|gy = B p • n|ay. In our case, where V is only bounded by S and extends 
to infinity above S (z > 0), and for A p , A, B p , B vanishing at infinity, the condition refers 
only to the boundary S, i.e. 

B z| 5 = B p z| 5 . (10) 

It is not necessary to independently pose equation (jTUJ), however, as it stems directly from 
equations (5) and (7b). Take the normal (vertical) components of equation (5a) and (5b) on 
S to obtain 

V fe -(A p xz)| (S = B p -z| <s , V h . (Ax z)\ s = B-z\ s , (11) 

respectively, where denotes differentiation on the horizontal plane, i.e. the boundary S. 
Equation ( jTOi) immediately follows from combining equations (7b) and ( jTTj) . In addition, 
from equations (J4]) and (JTUJ) we obtain B c • z|s = 0, so B c is a closed, toroidal magnetic field 
on S. In fact, B c is purely toroidal in any cross-section S' of V. This has been concluded 
by Kusano et al. (1994) and Berger (1999) who argued that the net toroidal flux of B p and 
B should be the same along any cross-section of V. The potential field B p being purely 
poloidal, on the other hand, one expects B p • B c = 0. A construction of B p and B c by 
poloidal and toroidal components, respectively, can also be found in Berger (1985). 



2.2. Magnetic energy and helicity budgets 

From the general equation (j3J) it is clear that the total magnetic energy E = [1/ (8tt)] j v B 2 dV 
of a closed magnetic structure is simply the sum of the potential magnetic energy E p = 
[l/(87r)] J v B 2 dV and the nonpotential magnetic energy E c = [l/(87r)] J v B 2 dV stored in the 
configuration in the form of electric currents: 

E = E P + E C . (12) 

Our objective will be to derive a convenient expression for each of the terms in equation 
( I12p . In this section we provide general energy expressions enabled by the gauge invariant 
definitions of the vector potentials A, A p , and A c . Equations for the potential energy E p 
can be directly applied to solar magnetic field measurements. Applicable expressions for the 
total energy E of the magnetic structure are given in the following sections, where the LFF 
approximation is adopted. 
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From the definition of A p , equation (5a), the potential magnetic energy of the configu- 
ration is given by 

E p = ^- A p x B p ■ n da , (13) 
87r JdV 

where da is the surface element on dV. To reach equation ([TBI we have used Gauss's theorem 
and the current-free condition, V x B p = 0. Since A p and B p , together with A and B, all 
vanish at infinity, however, the above surface integral applies only to the lower boundary S. 
The potential energy E p becomes, therefore, 

E p = I B p x A p • z dS . (14) 

Similarly, from the definition of A, equation (5b), the total magnetic energy of the 
configuration is given by (see also Berger 1988) 

E = — [ B x A p • z dS + — I A-VxBdV . (15) 
8vr Js 8tt J v 

The total energy from equation (|T5|) naturally tends to the potential energy in case the 
magnetic field vector B tends to its current-free limit B p . Decomposing B into B p and 
B c , one may derive equation ffT2"j) . where the potential energy is given by equation (|14p and 
E c =[l/(87t)]f v B 2 c dV. 

In view of the definitions and conditions for A and A p , equations (5)-(7), on the other 
hand, the relative magnetic helicity of equation (CQ) can be written as (e.g., Berger 1999) 

H m = / A • BdV , (16) 
Jv 

Equation (|14p enables the calculation of the potential energy E p for a flux-balanced 
magnetic configuration B, regardless of whether the magnetic field vector B is fully known 
on S. What is needed is the boundary condition for the vertical field B z = B ■ z on S, 
that uniquely determines the potential magnetic field B p and its vector potential A p on 
the boundary. In particular, assuming B p = —Vip, where ip is a smooth scalar, V 2 ip = 0, 
Schmidt (1964) showed that 

., . 1 f f B z (r')dx'dy' 

^' Z) = ^J J V(r-rr + *» ' ( > 

where r = xx + yy, r' = x'x + y'y are vector positions on S, defined for a given Cartesian 
coordinate system centered on S. For the vector potential A p one similarly obtains (see also 
DeVore 2000) 

POO 

A p (r,z)=Vxz *jj(r,z')dz' . (18) 
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Although exact, equations (fPT|) and fjlSp are computationally extensive. Much faster alter- 
natives are provided by means of Fourier transforms. Alissandrakis (1981), in particular, 
showed that 

B p (r) = F-\ j^== b u , v )k + F-\ ^= gy + B z z , (19) 
while Chae (2001) showed that 

A p (r) = ^(-j— - -2&«,„)* + F- x (^— ^6 Ui .)y , (20) 

where b U)V = YlaZi J2m=i B Zml exp[—i{ul + vm)] is the Fourier amplitude of B z , u = (2irl/L x ), 
v = (2nm/L y ), L x , L y are the linear dimensions of S, and F~ 1 (g) denotes the inverse Fourier 
transform of a function g. Albeit much faster, however, equations (|19p and (120]) assume 
periodic boundary conditions for B p and A p which contradicts the assumption of A and B 
vanishing at infinity. This problem is also well known. To mitigate the effects of the periodic 
boundary conditions assumed when Fourier transforms are used, one typically surrounds 
the initial flux concentration with a region of zero flux. In our calculations in §5.2 we have 
applied equations (fl9"l) and (120]) to zero-buffered magnetograms of solar active regions. 

From the above, E p can be readily calculated for flux-balanced photospheric or chro- 
mospheric (not necessarily vector) magnetograms of solar active regions. If a vector magne- 
togram is available, the vertical field B z on S is provided by rotating the measured magnetic 
field components to the local, heliographic, reference system (Gary & Hagyard 1990). Al- 
ternatively, the line-of-sight component can be used instead of B z provided that the studied 
active region is located sufficiently close to the center of the solar disk. This requirement 
typically minimizes the impact of viewing projection effects caused by the curvature of the 
solar surface. 

Unlike the potential magnetic energy E p , the total magnetic energy E, equation (TT5|) . 
cannot be calculated without additional assumptions or by using line-of-sight magnetograms. 
This is because A and V x B are generally unknown on and above the (photospheric or 
chromospheric) boundary S. 



3. The energy-helicity formula in the linear force-free approximation 

In the force-free approximation, (V x B) x B = 0, the total magnetic energy E of a 
structure extending in V is provided by the magnetic Virial theorem (Molodensky 1974; Aly 
1984) 

E = — [ \-B 2 R - (B ■ R)B1 ■ Ma , (21) 



where R is a vector position with arbitrary origin in V. For planar geometry, V extending to 
infinity and being bounded only by S at z — 0, and under the assumption that the magnetic 
field strength B vanishes with distance more rapidly than i? _3//2 , the Virial theorem reduces 
to its well-known form 



where r = xx + yy is a vector position with arbitrary origin on S. Equation (122]) has 
been applied to solar active regions (Metcalf, Leka, & Mickey 2005; Wheatland & Metcalf 
2006) assuming potential or force-free (not necessarily linear) magnetic fields. The explicit 
dependence of equation fl22|) on the coordinate system leads to inconsistencies if the employed 
magnetic field vector is not force-free (for a detailed discussion of problems related to the 
magnetic Virial theorem see Klimchuk, Canfield, & Rhoads 1992). Although well-known and 
particularly useful, the Virial theorem does not link the magnetic energy budgets with the 
relative magnetic helicity in a self-consistent way. For this reason, we will hereafter follow 
our alternative formulation for the potential magnetic energy (equation (fT4j) ) and the total 
magnetic energy in the LFF approximation. As shown in Figure [8] and explained in §5.2, 
the Virial theorem and our energy expressions give very similar results. 

Implementing the force-free approximation, V x B = aB, the total magnetic energy 
from equation (fT5l) gives 



and corresponds to the energy-helicity formula of Berger (1988). In case of the LFF approx- 
imation, where the force-free parameter a is constant in V, the dependence between E and 
the relative magnetic helicity H m becomes explicit. Substituting equation ffTBT) into equation 
( 123]) for constant a, we obtain 



Since the relative magnetic helicity depends entirely on the presence of electric currents so 
that H m = for B = B p , the first term in the rhs of equation (1241) must correspond to 
the magnetic energy that does not include the energy stored in electric currents for any 
nonzero a and H m . This ground-state energy can only be the potential energy E p , so the 
energy-helicity formula in the LFF approximation reads 




(22) 




(23) 




(24) 




m 



(25) 



A proof of equation (125]) is provided in Appendix A. 
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In a const ant- alpha magnetic structure the sense (sign) of the relative magnetic helicity 
H m is dictated by the chirality (sign) of the unique value of the force-free parameter a. 
Therefore, aH m > by definition in the LFF approximation, where a and H m are nonzero. 
If aH m = 0, then both a and H m are zero by definition. In this case, E = E p from 
equation (|25|) . By means of equation (fl2l) . moreover, the LFF approximation implies a 
linear dependence between the free magnetic energy E c and the relative magnetic helicity 
H m , namely, 

E c = ^-H m . (26) 

The monotonic dependence of E c on a and H m can be understood if one considers that both a 
nonzero relative helicity and a nonzero free energy depend on, and are directly proportional 
to, the existence of electric currents (a ^ 0). Of course, equation fl26l) is valid only for 
constant a because the sense of helicity is the same throughout the magnetic structure. For 
a non-constant a and in case of equal and opposite amounts of helicity being present in the 
structure, the net relative helicity H m becomes zero. This would give E c = for a ^ in 
equation ( |26i) . which is not true. 

Given that the potential energy is readily calculated (equation (fT4]) ). it is evident from 
equation (l25l) that knowledge of the relative magnetic helicity H m is sufficient to fully eval- 
uate the LFF energy-helicity formula. However, H m cannot be evaluated from the general 
equation (1161) . This is because the vector potential A is unknown in V, although the LFF 
magnetic field B can, in principle, be calculated everywhere in V. In the next section we 
derive convenient expressions for A and H m in the LFF approximation. 



4. The relative magnetic helicity in the linear force-free approximation 

4.1. Volume-integral representation 

It is straightforward to obtain a volume-integral expression for the relative magnetic 
helicity H m from the energy-helicity formula, equation (I25p . using the definitions of the 
potential and the total magnetic energies: 

H m = - !{B 2 - B 2 p )dV . (27) 
a Jv 

This expression was used by Hagino & Sakurai (2004) who assumed A = on S. Despite 
its simplicity, however, equation (T27j) cannot be directly compared to the general equation 
f|T6|) because the form of A is not obvious. To make this conceptual step, we decompose B 
in equation fT271) into its potential and nonpotential components to find, after some algebra, 
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that 

H m = - f B c • BdV , (28) 

where the condition B p ■ B c = has been used. Equation (1281) is identical to equation (1T61) 
for A = (l/a)B c , or 

A = -(B-B P ) . (29) 

a 

One may verify that the definition of A in the LFF approximation, equation (125)1 . 
complies with all the conditions of equations (5)-(7). Moreover, it is clear that A|s 7^ 0, 
contrary to what Hagino & Sakurai (2004) assumed. This being said, A does not have 
a vertical component on S because B and B p share the same vertical component on S, 
equation (TTOl) . Adopting A ■ z\$ = 0, however, equations (9) will give 

A c xz| 5 = A c -z| 5 = . (30) 

From equations (130)) . then, A c |s = 0, rather than A\$ = 0. Since A c |,s = 0, A p and A 
coincide on S via equation (jSJ), namely 

A| 5 = A P | 5 . (31) 

Of course, equation (13T)) does not preclude V x A|s 7^ V x A p |,s. This is because these curls 
include vertical gradient terms (d/dz) and A 7^ A p above S. 

Equation (1251) for the relative magnetic helicity satisfies all the requirements of the LFF 
approximation. Moreover, it is more complete than the equation H m = (1/a) J v B 2 dV of 
Pevtsov, Canfield, & Metcalf (1995). Clearly, from the above expression lim a ^o \H m \ = 
E p \im a ^ (l/\a\), which tends to infinity while it should tend to zero. Given that B in 
equation (1281) corresponds to an LFF magnetic field, the integrand B c • B = B 2 — B ■ B p 
is known at any location in V if B z is known and flux-balanced on S. Therefore, the total 
relative magnetic helicity of a closed and flux-balanced magnetic structure with constant 
a in V can be calculated using equation (1251) . By extension, equation (1251) can be used to 
calculate the total relative magnetic helicity of an isolated solar active region for which the 
constant-alpha approximation is assumed valid and for which flux-balanced photospheric or 
chromospheric vector magnetic field measurements exist. A representative, unique value of 
the force-free parameter a can be calculated by an array of techniques (Leka & Skumanich 
1999; Leka 1999) with an alternative technique described in §5.1. 

In theory, one may estimate the relative magnetic helicity from the volume integral of 
equation (1251) . In practice, however, it would be both risky and computationally extensive to 
apply equation (1251) to extrapolations of actual solar magnetic configurations. This is because 
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(i) despite the treatment of periodic effects on the horizontal plane, Fourier transforms in 
the LFF approximation can still cause periodic effects in the vertical direction and at large 
heights above the photosphere (i.e., Alissandrakis 1981; Gary 1989), and (ii) one does not 
know a priori the height above the photosphere where integration should stop. There are 
some partial remedies for both of the above problems: the use of a much larger computation 
volume than that required to contain the magnetic structure may limit periodic effects within 
V, while the maximum integration height can be either equal to the linear size of the surface 
S or determined by the contribution to the relative magnetic helicity H m . If integration 
above a certain height makes insignificant contributions to H m then integration stops at this 
height. In any case, calculating H m from equation fl28l is a very time-consuming task. For 
this reason, we derive in the next section a first-principles surface integral for H m in the LFF 
approximation. 



4.2. Surface-integral representation based on Fourier transforms 

Regardless of the force-free approximation, the total magnetic energy E of a closed 
magnetic structure extending into V and rooted in S consists of the potential magnetic 
energy E p of the structure and the nonpotential (free) magnetic energy E c due to electric 
currents (equation ffT2"j) ). Expressing E c in terms of E p , one may write 

E = (l + f)E p , (32) 

where E c = fE p and / is generally a positive and dimensionless variable. The constant- 
alpha approximation readily provides a condition for /, namely lim^o / — 0. In addition, 
/ must be a function of a and increasing |a| should increase / monotonically giving rise to a 
symmetric profile of f{a) with respect to a = 0, i.e. /(|a|) = f(—\a\). The form of / can be 
derived analytically in the LFF approximation if one uses the formulation of Berger (1985). 
The details of the derivation are given in Appendix B. Here we provide two expressions for 
the variable /. The first is the exact analytical formula, while the second is a linearized, 
with respect to a 2 , version f] of it, useful to keep the free energy and relative helicity finite 
when \a\d — > (2it/L), where d is the elementary size on the boundary S and L is the linear 
size of the magnetic structure on S. For an observed magnetogram, d corresponds to the 
linear size of a pixel expressed in physical units. In particular, 

/ = T , and fi = Tid^o? , (33) 

where 

L!=lZ J m=ir« Il ^l ( M 2 +B 2 i )l/2(„2 +1 ,2 i _ a 2 (i 2)l/2 

•r~ — [75 i j (34a) 

1=1 l^m=l ( M 2 + „2Jl/2 
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and 



1 2^1=1 Z^m=l ( u 2 +1) 2J3/2 





(34b) 



2^1=1 Z^m=l («2 +t ,2 j )i/2 



respectively. In equations (34), & Uj) „ m is the Fourier amplitude of the vertical magnetic field 
B z for the harmonic (ui,v m ) in a Fourier space with dimensions n x , n y . The linearization // 
implies a minimum value for / that results in the estimation of a minimum free magnetic 
energy E c and relative magnetic helicity \H m \ in the LFF approximation. The underestima- 
tion of E c and H m is negligible for small |a| and increases as a increases (see §4.3 below and 
Figured]). This effect is explained in detail in Appendix B and has been realized by several 
previous works when the linearized relative helicity expression of equation ([2]) was derived 
(Green et al. 2002; Demoulin et al. 2002; Demoulin 2006). The infinite energy and helicity 
for (|a|/c?) — > (2tt/L) is a well-known problem of the LFF magnetic fields that are not fully 
described by the boundary condition on S in this case (e.g., Alissandrakis 1981). 

The quadratic dependence of / on a in both the exact and the linearized case guarantees 
a symmetric profile of f{a) and a vanishing / for a — > 0. This dependence has also been 
demonstrated graphically by Sakurai (1981) in analytical force-free fields. 

From equations (I3"2"|) and (I3"3"j) we can now parameterize all terms of the energy-helicity 
formula (equation f|23]) ) with respect to one of these terms. As the free parameter we choose 
the potential magnetic energy E p for which the general expression of equation (fi~4")) exists. 
Then, the exact and the linearized surface-integral expressions for the total magnetic energy 
in a const ant- alpha magnetic structure read 





and 




(35b) 



respectively. For the free magnetic energy we obtain 



E c = J-E, 



'P 5 




and 



E c = Tid 2 c?E, 



'p i 



(36b) 



respectively, while for the relative magnetic helicity we find 

H m FEp , 

a 



(37a) 



and 



H, 



711 



%TiT i d 2 aE 1 



'p 1 



(37b) 



respectively. Notice that the exact formula for the relative helicity, equation (37a), still 
yields \H m \ — > for \a\ — > because T oc a 2 tends to zero faster than a. Equation (37b) 
gives values that are a factor of four smaller than those of the linearized relative helicity 
of Demoulin et al. (2002) and Green et al. (2002) (equation ([2])). This can be seen from 
equations (37b) and (1B3I) . and by setting d = l/(27r), which corresponds to Berger's (1985) 
unit length assuming a computational box with linear size L equal to unity. Part of the 
discrepancy has been corrected by Demoulin (2006) who admits that the original expression 
of equation ([2]) was a factor of two too high. In Appendix B we show that the linearization 
introduces an additional (l/2)-factor in equation (34b). 

Concluding, equations (35) - (37) offer convenient surface-integral representations of the 
relative magnetic helicity, as well as of the total and the nonpotential magnetic energies in 
the LFF approximation. Equations ( 1T4"1) and (35)-(37) will be evaluated in §§4.3 and 5 for 
semi-analytical and observed magnetic configurations, respectively. 

From equations (35) - (37) we can calculate the uncertainties SE, 5E C , and SH m of the 
total magnetic energy, the free magnetic energy, and the total magnetic helicity, respectively. 
Uncertainties of the potential energy, equation (1141) . stem from the uncertainties SB Z of the 
normal (vertical) magnetic field component B z . Although the values of 5B Z are generally 
known for a given magnetogram, it is difficult to propagate them into the potential energy 
because of the extrapolations required to infer the potential magnetic field and its vector 
potential. In our case the extrapolations are performed using Fast Fourier transforms. For 
this reason we will ignore the uncertainties SE P of the potential energy, although we expect 
that these uncertainties should not be very significant, given that the vertical magnetic field 
component is the least uncertain measured component of the magnetic field vector, especially 
for active regions located close to the center of the solar disk. For the same reason we will also 
ignore the uncertainties ST\ of T[. Excluding 5E P and STi, the only source of uncertainties 
is the uncertainties Sa in the inference of the force-free parameter a. These uncertainties 
give rise to a nonzero ST in the value of T (equation (34a)). From equations (35) - (37), 
then, we obtain the following uncertainty expressions: 

For the exact and linearized total magnetic energy, 



and 



SE 


E C ST 


(38a) 


~E 


~ ~E~T ' 


5E 
~E 


^E c Sa 
~E\a\ ' 


(38b) 



respectively. For the free magnetic energy, 

SE C ST 
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and 

# = 2^ , (39b) 
E c \a\ 

respectively. For the relative magnetic helicity, 



I -"ml V Ol J- 

and 

= ^ . ( 40b ) 
|-n m | M 

respectively. The "<" symbol in equation (40a) is due to the fact that T and a are 
interrelated. Given that it is also difficult to propagate the uncertainties of |a| into 5 J 7 , we 
will hereafter use the linearized expressions of the uncertainties, equations (38b) - (40b). 



4.3. Comparison between the volume- and the surface-integral expressions for 

the relative magnetic helicity 

At this point we have derived two types of expressions for the relative magnetic helicity in 
the LFF approximation, namely the volume integral of equation (1281) and the surface integrals 
of equations (37). To ensure consistency, these expressions must provide nearly identical 
results for small values \a\ of the force- free parameter, while the linearized expression of 
equation (37b) must provide a lower limit of the relative magnetic helicity as \a\ increases. 
To avoid errors due to observational uncertainties and to make a safer evaluation of the 
volume integral of equation (1251) we use semi-analytical models of magnetic structures, rather 
than observed solar magnetograms. For a simple representation of the twist present in the 
magnetic configurations we use dipolar magnetic field models. For a given dipole with 
footpoint separation L sep , we define the dimensionless quantity N = aL sep . This quantity is 
generally a dimensionless measure of a. In the particular case of field lines winding about an 
axis (not necessarily assumed here), N is a measure of the total end-to-end number of turns 
of the dipole. Our dipoles are characterized according to their iV-values. We first create 
the analytical distribution for the vertical magnetic field B z normal to the horizontal plane 
S and then we apply LFF extrapolations in the volume V using the same ^-distribution 
as boundary condition and assuming different a-values stemming from different iV-values in 
each extrapolation. Extrapolations are performed using the Fast Fourier transform method 
of Alissandrakis (1981). For this test we use positive a- values which results in right-handed 
helicities. Using equal and opposite a- values would only change the sense of twist and hence 
the sign, but not the magnitude, of the calculated magnetic helicity. The magnetic energy 
budgets, equations (35) - (36), are insensitive to the sign of a. 
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Our model dipoles have a fixed footpoint separation length L sep = 100 (in arbitrary 
units) and are embedded in a boundary plane S with linear dimensions L x = L y = 200, 
assuming an elementary length d — 1 as the unit length. The separation length L sep repre- 
sents the distance between the positive- and the negative-polarity centers of the dipole on 
S. We use an array of N- values, where N G [0.02,2], each of which determines a different 
a = (N/L sep ). To make sure that equations (35) - (37) do not depend on the details of a 
particular model, we use three different models of B z on S. For each model, the two polarity 
centers are placed at vector positions r x and r 2 , respectively, on S, such that |r x — r 2 | = L sep . 
The number of harmonics used for the Fourier-transform calculation of T and Ti, equations 
(37a) and (37b), respectively, is kept fixed in all cases and is equal to n x = n y = 256. We 
use the following models: 

(1) A Gold-Hoyle dipole solution (Gold & Hoyle 1960), i.e., 

*W = H + «(,'- ,,)» - ! ' <41) 

where r is the vector position of a given location on S and B , q are positive constants. 
In this test we have used a fixed B = 10 3 and an array of q- values, q G [0.5, 10]. Each 
g-value has been applied to the full array of iV-values. 

(2) A solenoidal dipole solution (Sakurai & Uchida 1977), i.e., 

1671 tt ^16(1 +^2 + 7^ 16(1 ~ Pi) +nL °<* 

(42) 

where Si = 1, s 2 = —1, q is a positive constant, pi = (pi/q); % = {1,2}, are the 
normalized, with respect to q, distances of a location r on S from i"! and r 2 (pi = |r — r*ij ; 
i = {1,2}), L sep = L sep /q is the normalized, with respect to q, separation length, 
kf = 64pi/[16(l +pi) 2 + Tc 2 L 2 sep ], and h(ki), hih) are the complete elliptic integrals of 
the first and second kind, respectively, i.e., 

rw/2 jq Ptt/2 , 

h(k t )= and I 2 (k t )= Jl-kfsinme . (43) 

Jo a/1 — kfsm 2 9 Jo v 

Here we have used a fixed B = 10 3 and an array of q- values, q G [0.02, 0.3], each value 
of which has been applied to the full array of N- values. 

(3) A submerged poles dipole solution (Longcope 2005 and references therein), i.e., 

Bz{r) = B ° q{ [(r - + g2]3/2 " [(r _ ^ + g2 ]3/2 > ' ( 44 ) 
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where Bq and q are positive constants. The constant q, in particular, represents the 
depth below S in which the two magnetic monopoles are placed. The depth of each 
monopole can, in principle, be different than that of the other (s), but here we use a 
fixed depth, as well as a fixed magnetic field strength B for each monopole, to create 
a flux-balanced magnetic configuration on S. Here we use Bq = 10 3 and an array 
of depths q, q e [0.5,10], each value of which has been applied to the full array of 
iV-values. 

Comparing the volume-integral expression, equation (12811 . with the surface-integral ex- 
pressions, equations (37), for the relative magnetic helicity H m gives the expected results 
for all q- and iV-values. Three of these results, one for each model, are given in Figure 
[Ik. The q- values for each model in Figure [Th^ were selected with the sole purpose of giving 
rise to well-separated helicity values, for convenience in the visual comparison. These selec- 
tions are q = 0.6, 0.2, and 10 for the Gold-Hoyle (GH), the Sakurai-Uchida (SU) and the 
submerged-poles (SP) model, respectively. Figure [Th, gives rise to the following conclusions: 

(1) All helicity expressions give very similar results for a given model, which suggests that 
the LFF equations (35) - (37) are model-independent. 

(2) For iV = aL sep — > 0, all expressions give H m — »■ 0. Therefore, H m from equations 
(J2"8"j) and (37) corresponds to the gauge-invariant relative magnetic helicity discussed 
in §§4.1 and 4.2. As a increases for a fixed L sep , Figure [Jji shows the quadratic increase 
of the magnetic helicity in the dipoles, for a fixed boundary condition B z on S. 

(3) Clearly, all expressions for H m give nearly identical results for small a. As a increases, 
the linearized surface-integral expression (equation (37b); dotted curves) consistently 
provides a lower H m (the exact surface-integral expression of equation (40a) is repre- 
sented by solid curves and rectangles), as expected. The volume- integral expression 
for H m (equation (12"8"|) ; dashed curves and triangles) gives slightly higher values than 
both surface-integral expressions. 

Tests with N > 2 (not shown in Figure [T]) reveal that the volume-integral H m increases 
exponentially after some maximum iV-value, while the surface-integral expressions continue 
to increase quadraticalljIE This maximum iV-value is model-dependent and, in case of the 
SP model, it changes even with varying model parameters. This inability to predict the 



2 Of course, at the vicinity of \a\d ~ (2ir/L), which in our parameter selection corresponds to N — 7T, the 
exact surface-integral H m increases abruptly to become infinite for \a\d — (2n/L) 
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behavior of the volume-integral H m for large N, combined with spurious results of the 
Fourier-transform extrapolations in these cases, enhances one's impression that the volume- 
integral helicity is less reliable than the surface- integral expressions and that, among other 
problems, it is susceptible to artifacts incurred by Fourier-transform extrapolations at large 
heights above the boundary S. 

Figure [lb shows the underestimation factor caused by the use of the linearized surface- 
integral expression of equation (37b) for the three cases depicted in Figure dh- In particular, 
the dashed curves show the ratio between the volume-integral and the linearized surface- 
integral helicities, while the solid curves show the ratio between the exact surface-integral 
and the linearized surface-integral helicities. Evidently, the underestimation factor is nearly 
model-independent when the two surface-integral expressions are compared, while the situ- 
ation is less predictable when the volume-integral and the linearized surface-integral expres- 
sions are compared. From the comparison between the surface-integral helicity expressions, 
one sees that, even for large iV-values, underestimation does not exceed a factor of ~ 1.85. 
For dipolar solar active regions for which the LFF approximation is assumed and with typical 
values of a ~ 1CT 2 Mm and L sep ~ 100 Mm, we obtain N ~ 1. The expected underesti- 
mation factor for this case is < 1.1, which is very modest compared to the errors expected 
from other assumptions, and especially the use of the constant-alpha approximation itself. 

In summary, the results shown in Figure [1] demonstrate that the surface-integral expres- 
sions of equations (37) lead to reliable estimates of the total relative magnetic helicity in a 
constant-alpha magnetic structure. By extension, the surface integrals of equations (35) and 
(36) provide reliable estimates of the total magnetic energy and the free magnetic energy, 
respectively, in the structure. 

5. Application to observed solar active region magnetic fields 

5.1. Data selection and determination of basic parameters 

In this section we apply the results of the previous analysis to vector magnetograms of 
solar active regions. In particular, we calculate the LFF magnetic energy and helicity bud- 
gets (equations (TH|) and (35) - (37)) using photospheric vector magnetogram data obtained 
by the Imaging Vector Magnetogram (IVM; Mickey et al. 1996; LaBonte, Mickey, & Leka 
1999) of the University of Hawaii's Mees Solar Observatory. IVM's photospheric magnetog- 
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raphjo consists of recording the complete Stokes vector at each of 30 spectral points through 
the Fe I 6302.5 A photospheric spectral line. The line-of-sight magnetic field components are 
obtained via the inversion code of Landolfi & degl'Innocenti (1982) that includes LTE radia- 
tive transfer, magneto-optic effects, and the filling factor of the unresolved flux tubes. Our 
equations are applicable to the heliographic magnetic field components on the heliographic 
plane. The required coordinate transformation has been carried out following the analysis 
of Gary & Hagyard (1990). 

To test our derivations for energy and helicity, we study two solar active regions (ARs): 
a small, short-lived emerging flux region (NOAA AR 8844) and a persistent, large, and 
complex AR that exhibited significant eruptive activity (NOAA AR 9165). For both ARs, 
the IVM recorded a series of vector magnetograms over a period of a few hours that can be 
used to follow the temporal evolution of the energy and helicity budgets in the regions. 

Before applying any analysis that employs the transverse field of a vector magnetogram, 
one must first resolve the intrinsic azimuthal ambiguity of 180° in the orientation of the 
transverse magnetic field component. Azimuth disambiguation of the employed IVM mag- 
netograms was performed by means of the nonpotential magnetic field calculation (NPFC) 
method of Georgoulis (2005a) - see also Metcalf et al. (2006) for a comparative evaluation 
of the method with respect to other disambiguation methods. Figures [2^ and [2Jd depict two 
disambiguated vector magnetograms of NOAA ARs 8844 and 9165, respectively. Only part 
of the IVM field of view is shown in both images, to exemplify the magnetic structure of the 
two ARs. Shown are the heliographic magnetic field components on the heliographic plane. 
The relative isolation of the two ARs on the solar disk at the time of the IVM observations 
(2000 January 25 and September 15 for ARs 8844 and 9165, respectively), as well as the 
ARs' very different records of eruptive activity prompted us to use these data in this first 
test of our LFF energy / helicity calculations. 

The timeseries of the magnetic flux $ during the IVM observing interval for both ARs are 
shown in Figure [31 For both cases, we notice that the IVM field of view encloses fairly well- 
balanced magnetic flux distributions. NOAA AR 8844 is more flux-balanced than NOAA 
AR 9165, with an imbalance always kept below 5%. The maximum imbalance of NOAA 
AR 9165 is around 10%. Our derivations require flux-balanced magnetic structures and 
the above slight imbalances are not expected to significantly impact our results. The first 
noticeable difference between the two ARs is in their respective amounts of magnetic flux: 
on average, the magnetic flux in NOAA AR 9165 ($ ~ 17.1 x 10 21 Mx) is a factor of ~ 3.4 



3 Recently, the IVM focused on the chromospheric magnetically sensitive line Na I (5896 A). These ob- 
servations have started providing chromospheric vector magnetograms of solar active regions. 
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larger than the flux in NOAA AR 8844 ($ ~ 5.1 x 10 21 Mx). One might also notice a very 
slight increasing trend in the evolution of $ in NOAA AR 8844 (from ~ 4.9 x 10 21 Mx to 
~ 5.3 x 10 21 Mx) within the 2 hr of the IVM observations, implying that the magnetic 
structure is growing. This is typical of emerging flux regions. 

After disambiguation, we need to calculate a unique force-free parameter a for each 
magnetogram. To do so, we calculate the slope in the scatter plot between the vertical curl 
(V x B) z of the magnetic field B and the vertical field B z for strong-field locations of the 
magnetograms. By strong-field locations we mean locations with magnetic field components 
exceeding the la threshold, where we have taken lcr to correspond to a vertical magnetic 
field of 100 G and a horizontal magnetic field of 200 G, typical of the IVM. Because the 
LFF approximation is a gross simplification of the photospheric active-region magnetic fields, 
however, the uncertainty in the value of the slope is often larger than the slope itself due to 
the substantial scatter in the pairs of [(V x B) 2 , Bj-values. To restrict the uncertainty in 
the calculation of a we have developed the following procedure: we obtain several a-values 
from the slope of the scatter plot, each calculated using a different significance threshold, as 
shown in Figure HI Let a^ be the value of a for a given significance threshold = ka; k > 1 
and be the respective unsigned magnetic flux used in the calculation. Then, the adopted 
unique value of a and its uncertainty 5a are obtained by the flux-weighted averages 

a = E * Q f k and fa = ^> ' ^ , (45) 

respectively. An example of this calculation is shown in Figure HI The flux- weighted average 
a is indicated by the solid line and the surrounding shaded area indicates the extent of its 
uncertainty 5a. The above process provides a maximum-likelihood a-value with a reasonable 
uncertainly and is repeated for every vector magnetogram of the timeseries to obtain the 
respective timeseries for a. Using other methods to calculate a (see, e.g., Leka & Skumanich 
1999 and Leka 1999) we verified that the timeseries of a obtained by equations (I45p are more 
smooth (less spiky) and with smaller uncertainties for each a-value, than the a-timeseries 
stemming from the other methods. 



5.2. Magnetic energy and helicity calculations 

The timeseries of the force- free parameter a for both tested ARs are shown in Figure [5J 
As we discussed in §5.1, the different a-values are generally consistent with each other, giving 
rise to fairly well-defined averages in both cases. The overall twist for NOAA AR 8844 is 
right-handed (a > 0), while for NOAA AR 9165 it is left-handed (a < 0). For the latter AR, 
in particular, a appears to decrease, in absolute value, in the course of time. Coincidentally, 
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the average absolute values |a| of a for both ARs are almost identical (a = 0.023±0.06 Mm 1 
and a = -0.024 ± 0.06 Mm 4 for NOAA ARs 8844 and 9165, respectively). We note in 
passing that the value of a for NOAA AR 8844 is in excellent agreement with the value 
of 0.022 Mm -1 , calculated by Pariat et al. (2004). The latter a- value was inferred by 
combining the best LFF match of the active-region corona using simultaneous EUV images 
from TRACE with a best LFF fit of the observed horizontal magnetic field. The magnetic 
field vector in Pariat et al. (2004) was acquired by the high-resolution vector magnetograph 
onboard the balloon-borne Flare Genesis Experiment (FGE; Bernasconi et al. 2001). 

Although the two studied ARs happen to have almost the same a-values, albeit with 
different signs, the much larger magnetic flux carried by the eruptive NOAA AR 9165 is 
expected to lead to much larger energy / helicity budgets than the respective budgets of the 
noneruptive NOAA AR 8844. The relative magnetic helicity and the respective magnetic 
energies for ARs 8844 and 9165 are plotted in Figures [6] and [7J respectively. There we 
show both the linearized (equations (35b)-(37b); red curves) and the exact (equations (35a)- 
(37a); blue curves) surface- integral expressions for energy and helicity. Before discussing 
and comparing individual values, we note that the linearized expressions generally provide 
slightly lower magnitudes of energy and helicity. This is clearly the case for NOAA AR 
8844 (Figure [6]), while in cases where the linearized values are larger than the exact values 
for NOAA AR 9165 (Figure [7]), the difference is within error bars. That the linearized 
expressions provide lower limits of energy and helicity was concluded from the analysis in 
Appendix B and verified using semi-analytical models in §4.3. In both cases of observed 
ARs, moreover, the timeseries of the exact values appear more spiky than the respective 
timeseries of the linearized values. Given that the linearized expressions lead to a smoother 
temporal evolution with slightly lower values than the exact expressions, equations (35b) - 
(37b) for the linearized energy and helicity budgets appear preferable compared to the exact 
expressions of equations (35a) - (37a). Besides being more convenient and well-behaved, 
the linearized expressions also have readily derivable uncertainties (equations (38b) - (40b)) 
based on the uncertainties 5a of a. This being said, one notices the close correspondence of 
the a- value timeseries of Figures [5] with the timeseries of the linearized helicity of Figures [6] 
and [7J Clearly, the success of the LFF energy/helicity estimations depends on the reliability 
of the inference of a. This is a key feature that one should keep when trying to generalize 
the LFF energy/helicity formulas into NLFF ones, valid for a variable a within the field of 
view. 

The average energy /helicity values from Figures [6] and [7J are summarized and compared 
in Tables [T] and [2J Table [T] shows the comparison between average magnetic fluxes, a- values, 
and helicities, while Table |2] focuses on the comparison between the various average energy 
budgets from the two ARs. It is quite useful that the average a-values are nearly identical 
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for the two ARs. We then notice that the eruptive NOAA AR 9165, with a factor of ~ 3.4 
larger magnetic flux than the noneruptive NOAA AR 8844, also has a potential and a total 
magnetic energy that are similarly (by a factor of ~ 3.2 — 3.5) larger than those of NOAA AR 
8844. The average relative magnetic helicity and free magnetic energy of the eruptive AR, 
however, are ~ 7.6—8.9 times larger than those of the noneruptive AR. Notably, the linearized 
expressions, that result in lower values and uncertainties, consistently give a higher factor of 
difference in both energy and helicity. The much wider difference between the free energy 
and the relative helicity between the two ARs suggests that a viable criterion (safer than 
simply evaluating the magnetic flux) for distinguishing between eruptive and noneruptive 
ARs may be the amount of free magnetic energy and helicity. This study should, of course, 
be applied to a large number of eruptive and noneruptive ARs for the results to obtain 
statistical significance. Moreover, if this treatment is generalized to the NLFF, rather than 
the LFF, approximation, it will be much more physically meaningful given the expected 
conditions in the low solar atmosphere. Forced photospheric fields (see, e.g., Georgoulis 
& LaBonte (2004)) should still lead to discrepancies stemming from the application of the 
force-free approximation. Another notable fact from Table |2] is the fractional free magnetic 
energy (E c /E), normalized by the total magnetic energy. For the noneruptive NOAA AR 
8844, the free energy is ~ 4.7% - 6% of the total magnetic energy. For the eruptive NOAA 
AR 9165, the free energy corresponds to ~ 12.3% - 13.1% of the total energy, which is a 
factor of ~ 2.2 — 2.6 higher. 

Notice the the above ratios of the free to the total energy are substantially lower than 
those calculated by Metcalf, Leka, & Mickey (2005), for NOAA AR 10486, on 2003 October 
29. Using the Virial theorem and implicitly assuming NLFF magnetic fields, these authors 
found that the free energy ranged from ~ 44% (a few hours before a major X10 flare) to 
~ 75 — 80% (in the course of, and shortly after, the flare), of the total energy. These ratios 
appear extraordinarily high, at least in view of eruption models that predict the eruption 
onset when the free energy exceeds 10 — 15% of the total energy (see, for example, DeVore 
& Antiochos 2005). Of course, NOAA AR 10486 was an extraordinary AR, which might 
account for its unusual behavior. 

Despite the large difference of energy and helicity budgets between the two ARs, notice 
that significant magnetic helicity is present even in the noneruptive NOAA AR 8844. Indeed, 
the average relative helicity of the AR is H m ~ (1.5 ± 0.4) x 10 42 Mi 2 , with the helicity 
budget of a typical CME estimated at ~ 2 x 10 42 Mx 2 (DeVore 2000). With a minor helicity 
increase, therefore, the AR should be capable of producing a typical CME before relaxing to 
the potential state. Interestingly, a faint halo CME occurred above the AR on 2000 January 
26 at ~ 12:00 UT and the AR started to decay less than 24 hr later, on 2000 January 27 
(Schmieder et al. 2004). As the AR was still growing during the IVM observations, it is 
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likely that its magnetic helicity was further increased by January 26. No significant flaring 
activity was associated to the CME. 

NOAA AR 9165, on the other hand, gave an eruptive M2 flare a few hours before the 
IVM observations on 2000 September 15, as well as two even stronger eruptive flares (M5.9 
and M3.3) on the next day. Its relative magnetic helicity, ~ (—13 ± 4) x 10 42 Mx 2 was 
enough to launch nearly seven typical CMEs. Perhaps not surprisingly, the AR survived for 
several more days and could clearly be followed until it crossed the western solar limb. 

Back in our analysis, the compromise brought by the LFF approximation is reflected 
on the uncertainties accompanying the estimations of the free magnetic energy in both ARs. 
Obviously, the free magnetic energy is a crucial parameter in assessing the eruptive potential 
of a given AR (see, e.g., Metcalf, Leka, & Mickey 2005). With average linearized free energies 
of (0.15±0. 1) x 10 32 erg and (1.33±0.6) x 10 32 erg for NOAA ARs 8844 and 9165, respectively, 
the lowest uncertainties are estimated at ~ 67% and ~ 45%, respectively. A generalization 
allowing NLFF fields will hopefully restrict these uncertainties. 

Finally, in Figure [S] we compare our estimated potential and total magnetic energies, 
equations (fT4"l) and (35) - (36), with those obtained by the calculation of the Virial theorem, 
equation fl22|) . for both ARs. To implement the Virial theorem, we perform a current-free 
and a LFF extrapolation of each IVM magnetogram, the latter using the inferred maximum- 
likelihood a-value for this magnetogram. Estimates of the Virial theorem are represented by 
dashed curves and triangles. Solid curves and rectangles refer to the linearized expressions, 
while dotted curves refer to the exact expressions. For a convenient comparison, the scaling 
for the potential energy (blue curves) is different than the scaling for the total energy (red 
curves). From the plots in Figure [BJ we first notice that our potential-energy expression, 
equation (114)) . gives almost identical results with the Virial theorem for both ARs. The 
average fractional differences \E V — EL„. . J /(-En + E v .„. . ,,) are ~ 0.7% and ~ 0.2% for 
ARs 8844 and 9165, respectively. The difference is larger for the total energies. On average, 
the fractional difference is ~ 1.6% (~ 1.1%) for the linearized (exact) expressions in NOAA 
AR 8844. For NOAA AR 9165, the average fractional difference is ~ 2.7% (~ 3.2%) for the 
linearized (exact) expressions. These differences are small and generally within the uncer- 
tainties in the calculation of energies. In addition, the Virial theorem provides consistently 
slightly higher total energies for NOAA AR 8844, while it consistently provides slightly lower 
total energies for NOAA AR 9165. This mixed behavior, as well as the source of the slight 
discrepancy in total energies, are unclear. One possible reason may be the application of 
Fourier transforms, and hence an implicit assumption of periodic boundary conditions, in 
analytical expressions where fields are required to vanish at infinity. Why this does not have 
an impact in the calculation of the potential energy is also unclear. Nevertheless, the small 
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discrepancies prompt us to conclude that our energy expressions are consistent with the 
Virial theorem. The reasons why we have derived and used them instead of the latter are 
that (i) they provide a self-consistent description of the energy and helicity budgets, and (ii) 
they are physically intuitive, derived from first principles, and, hopefully, capable of being 
generalized for NLFF magnetic fields. 

6. Summary and discussion 

The reliable calculation of the magnetic energy and helicity budgets in the active-region 
solar corona is an essential step toward the quantitative understanding of solar eruptions and 
has profound space-weather applications. Our goal is to derive a practical set of equations 
that are applicable to solar vector magnetograms and can evaluate the magnetic energy and 
relative magnetic helicity budgets in a physically intuitive, self-consistent manner. Here we 
provide expressions for the magnetic energy and relative helicity budgets in case of a constant- 
alpha, flux-balanced, magnetic structure, thus implementing the LFF approximation. These 
equations are to be generalized into magnetic structures with non-constant alpha values, thus 
implementing the NLFF approximation. This objective will be pursued in a later study. 

To perform our LFF analysis we separately derive each of the terms present in the 
energy-helicity formula of Berger (1988), namely the total magnetic energy, the potential 
magnetic energy, and the relative magnetic helicity related to the free magnetic energy, to- 
gether with their uncertainties. Our analysis unifies numerous expressions for the relative 
helicity and links several virtually unconnected studies into a self-consistent energy-helicity 
description that is practical enough to be applied to vector magnetograms of solar active re- 
gions. For the ground-state, potential, magnetic energy we provide a general surface-integral 
expression, equation (TTJJ. This expression gives results practically identical to those of the 
magnetic Virial theorem. The potential magnetic energy is then used as a free parameter 
to explicitly determine the total and free magnetic energy, as well as the relative magnetic 
helicity. The variable relating the potential energy to the free energy and the relative helicity 
has been calculated in two ways - an exact and a linearized one - by using and extending 
the analysis of Berger (1985). As a result, the magnetic energy and helicity budgets are 
calculated self-consistently as surface integrals, equations (35) - (37). This development 
reduces significantly the required computations. Reliability and computational speed are 
essential elements of a future real-time or near real-time calculation of the magnetic energy 
and helicity budgets in solar active regions. 

To test our derivations we used three different types of semi-analytical LFF magnetic 
dipoles (§4.3). The conventional volume-integral expression for the relative magnetic helic- 
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ity was compared with our exact and linearized surface-integral expressions. The convincing 
match between the volume- and surface-integral helicity expressions for all models implies 
that our formulations are model-independent. Moreover, as expected from the analysis, the 
linearized surface-integral expression of the relative magnetic helicity, equation (37b), con- 
sistently provides a lower limit of the helicity present in the magnetic structures, with this 
behavior being more pronounced for large alpha values. For smaller alpha, all helicity expres- 
sions give nearly identical results. By extension, the linearized surface-integral expressions 
provide reasonable lower limits of the total and free magnetic energy, equations (35b) and 
(36b), respectively. Given also the convenient calculation of uncertainties in the linearized 
case, equations (38b) - (40b), we conclude that linearization is preferable over using the 
exact formulas, where alpha can resonate with the value of (2n/L) and hence lead to infinite 
energies and helicities, as is well-known for LFF magnetic structures. 

Two series of solar vector magnetograms, one for an eruptive and another for a nonerup- 
tive active region, were thereafter subjected to our analysis (§5). Both the exact and lin- 
earized expressions for the energy and helicity were used. We found that the exact expres- 
sions tend to give more spiky temporal evolutions and hence larger uncertainties in temporal 
averages, which provides an additional reason for preferring the linearized energy/helicity 
expressions in observations. Both ARs happened to exhibit nearly the same absolute alpha 
value. The eruptive active region, however, included several times more magnetic flux than 
the noneruptive AR. This effect was greatly amplified when the free magnetic energies and 
relative magnetic helicities of the two active regions were calculated and compared. This 
leads us to the conclusion that comparing the free energies and helicities might be a safe 
(safer than simply calculating the total magnetic flux) way of distinguishing between erup- 
tive and noneruptive active regional. The crucial point, however, is the reliable calculation of 
free energies and helicities. The LFF approximation is certainly not very reliable, as can be 
seen from the large error bars accompanying our free energy and helicity estimates (~ 45% 
- 70%). To reach sound conclusions, the analysis involving and comparing free magnetic 
energies / relative magnetic helicities must be applied to statistically significant samples of 
active regions, and ideally by utilizing the NLFF approximation. 

The force-free approximation is a prerequisite for our analysis because it is a very dif- 
ficult, if not intractable, problem to perform non-force-free calculations of the nonpotential 
magnetic energy and the magnetic helicity in active regions. The only hope for non-force-free 
energy and helicity calculations emerges from data-driven three-dimensional magnetohydro- 
dynamical (3D MHD) simulations of the active-region corona (e.g. Abbett 2003; Roussev et 



4 For an alternative criterion, based on the magnetic connectivity in solar active regions, see Georgoulis 
& Rust (2007). 
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al. 2004) which, however, require immensely time-consuming calculations. 3D MHD models 
are certainly capable of advancing our physical understanding of solar eruptions but, because 
of their intense computations, they cannot contribute to a real-time, or near real-time, space 
weather forecasting capability. Alternatively, non-force-free energy and helicity estimates 
can be obtained if an active region is continuously observed from its formation and there- 
after. In this case, total energies and helicities can be calculated by temporally integrating 
the Poynting flux and magnetic helicity injection rate, respectively. If the birth of an active 
region is not observed, then the initial energy and helicity can only be assumed. In any 
case, both the Poynting flux and the helicity injection rate require the flow velocity of the 
magnetized plasma on the boundary of the magnetic field measurements. Inferring a reliable 
flow velocity is a completely independent, as well as highly nontrivial, problem (for a review, 
see Welsch et al. 2007). 

Our force-free equations are physically better suited to apply to chromospheric, rather 
than photospheric, vector magnetograms. It has yet to be established whether the NLFF 
approximation holds for the active-region chromosphere (see Metcalf et al. [1995] in con- 
junction with Socas-Navarro [2005]) but it is almost certainly more valid there than in the 
photosphere. The first high-quality chromospheric vector magnetograms have already been 
obtained (the above authors as well as Leka & Metcalf 2003; Metcalf, Leka, & Mickey 2005; 
Wheatland & Metcalf 2006) but a routine acquisition and reduction of such data may still 
be a task for the future. In brief, force-free equations may be applied to photospheric vector 
magnetograms as a zero-order (LFF) or first-order (NLFF) approximation, but one expects 
larger uncertainties in the values of energy and helicity budgets, than when chromospheric 
vector magnetograms are used. 

Concluding, we emphasize that the present analysis cannot fully uncover the importance 
of magnetic helicity in solar eruptions. Here we only show two examples that appear to point 
to this direction but an answer would require large numbers of active regions and NLFF 
energy/helicity equations, as already said. Our objective here was to calculate the relative 
magnetic helicity in active regions as an integral part of the energetics and complexity of the 
studied magnetic structures. It would be an important leap forward if it was convincingly 
shown that flare- and CME-prolific active regions exhibit significant quantitative differences 
in their free magnetic energy and/or total relative helicity (large free magnetic energy does 
not necessarily imply a large total relative magnetic helicity because roughly equal and 
opposite amounts of helicity may be simultaneously present - see Phillips, MacNeice, & 
Antiochos [2005]) compared to quiescent ARs. Intriguing clues to this direction stem from 
the study of the structural magnetic complexity in solar active regions (e.g. Georgoulis 
2005b; Abramenko 2005) or the calculation of the free magnetic energy in active regions 
with exceptional flare and CME records (Metcalf, Leka, & Mickey 2005) but the role of 
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helicity is yet to be uncovered. Some pieces of evidence suggesting the importance of helicity 
in solar eruptions stem from the frequent presence of sigmoids in eruptive active regions 
(Rust & Kumar 1996; Canfield, Hudson, & McKenzie 1999), apparently due to significant 
amounts of helicity with a prevailing sign, the presence of large and highly variable alpha 
values in eruptive active regions (Nindos & Andrews 2004), and the statistical correlation 
between large helicity injection rates and X-class flares/CMEs (LaBonte, Georgoulis, & Rust 
2007). Our forthcoming NLFF analysis will be well suited to address the role of helicity in 
solar eruptions and we intend to carry out this study in the future. 

This work is dedicated to the memory of its co-author, Barry J. LaBonte. Barry is 
remembered as a deeply knowledgeable, distinguished colleague and an inspiring mentor. I 
am grateful to D. M. Rust for our continuous interaction on magnetic helicity in the Sun and 
for a critical reading of the manuscript. I also thank A. Nindos and S. Regnier for clarifying 
discussions on helicity issues and an anonymous referee whose numerous critical comments 
and suggestions resulted in substantial improvements in the paper. Partial support for this 
work has been received by NASA Grants NAG5-13504 and NNG05-GM47G. 

A. Equivalence of equations (124ft and ( 1251) for the energy-helicity formula in 



To show that equations (1241) and (1251) are equivalent in the LFF approximation, it is 
sufficient to show that the potential energy E p is given by 



for any LFF magnetic field B ^ B p . 

We first decompose B in equation (lAlj) into its potential (poloidal) and nonpotential 
(toroidal) components, B p and B c . Then, equation (lAlj) becomes 



The first integral of equation (1A2I) is already the potential energy as shown in equation (fl4|) . 
To prove equation (1A1I) . therefore, it is sufficient to show that 



the linear force-free approximation 




(Al) 




(A2) 




(A3) 
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From Bp • B c = 0, we construct the volume integral J v B p • B c <iV = 0. Substituting the 
definition of A p from equation (5a) into this volume integral, we find after some analysis 
that 



A p x B c ■ nda = - / A p ■ V x B c dV . (A4) 

dV JV 

Taking into account that (i) A p vanishes at infinity, and (ii) V x B c = V x B, because 
V x Bp = 0, equation (]A4j) further reduces to 



B c x A p • MS = - / A p • V x BdV . (A5) 



In the LFF approximation, however, V x B = ctB, with a constant, so equation flA5j) gives 



B c x A p • MS = -a / A p ■ BdV . (A6) 



iv 



Given the gauge-invariant definition of the relative magnetic helicity, however, it can be 
shown (Berger 1988; 1999) that 

' A p • BdV = . (A7) 



v 



Combining equations (IA6I) and (IA7I) we obtain equation flA3j) . Therefore, equation flAlj) is 
true and hence equations (T24"|) and (j25|) in §3 are equivalent as asserted. 



B. Derivation of the variable linking the potential and the total magnetic 

energy in the LFF approximation 

Here we will derive the form of the dimensionless variable / in equation fl32|) . This 
variable links the total and the free magnetic energies in a const ant- alpha magnetic structure. 
We will use and extend the analysis performed in Appendix All of Berger (1985). Assuming 
planar geometry, Berger (1985) utilized Chandrasekhar's (1956; 1961) decomposition of an 
arbitrary magnetic field vector into a poloidal and a toroidal components and, in view of the 
LFF approximation, he derived the total magnetic energy and the relative magnetic helicity. 

Following Berger (1985), the total magnetic energy of the structure is given by 

n x n v I t 2 I 

E = lJ2J2-f^ > ( B1 ) 

Z 1=1 m =l H ' m 

where b UuVm is the Fourier amplitude of the vertical magnetic field B z for the harmonic 
(ui,v m ) in a two-dimensional Fourier space with linear dimensions n x , n y . In addition, we 
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have kf m = uf + v ^ — a' 2 . The force-free parameter a' is expressed in inverse length units 
(i.e., 1/x, where x is the number of unit lengths required for a' = 1) and not in physical units. 
This is why it is represented by a', while the a used so far refers to the force-free parameter 
expressed in physical units. Typically, a' = ad, where d is the unit length expressed in 
physical units. Berger (1985) assumes periodic boundary conditions and a length unit of 
[L/(2tt)\, where L is the linear dimension of the magnetic structure on the boundary S. 
Moreover, ui = (2ttI/L) and v m = (2nm/L). Then, the direct and inverse Fourier transform 
of B z can be performed on S only so that one can write 

Hi n y 

B z {x,y)\ s = J2J2 K* m e l{uiX+VmV) ■ (B2) 

1=1 m=l 

The required boundary conditions for b Ul)Vm in order to have a real and finite magnetic energy 
and helicity is b Ul , Vm = for \/u 2 + v 2 < \a'\ (see also Alissandrakis 1981). 

Assuming that the magnetic structure does not include electric currents (a' = 0), then 
equation (IB1I) provides the potential magnetic energy of the structure, namely 



^ = IEE^> (B3) 

Z 1=1 m =l 



b 2 

ui ,v m l 



J 2 . 



where q 2 l m = u] + v\ 

From equation the variable / is given by the dimensionless ratio 

/ = ^. (B4) 



Substituting equations fIBlj) and flB3j) into equation (1B4I) we obtain 



152 \<r± 

„ i—iu z—/v I u,v\ kq / nr \ 
/ = |62 I > ( B5 ) 



^ til ^ 



where we have denoted Ym=i E^=i ^y J2 U J2 V f° r simplicity. The ratio of sums in equation 
(1B5[) depends on a' because of its dependence on k. This dependence can cause problems 
when I a' I — > (2tt/L) because k — > for I — m = 1 in this case and / becomes infinite. This 
problem is not new; that LFF fields sometimes give solutions that are not fully specified by 
the boundary condition and may include infinite energy has been explicitly acknowledged 
by Alissandrakis (1981), but also by Chiu & Hilton (1977), using a different analytical 
framework. Clearly, this is a caveat of the LFF approximation and restricts its applicability. 
To avoid infinite energy values when \a'\ — > (2n/L), equation (IB5I) can be linearized with 
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respect to a /2 , assuming small values of a'. We first write (q — k)/(kq) = (q 2 — kq)/(kq 2 ). 
Expanding q 2 — kq in a MacLaurin series, one finds q 2 — kq ~ (q 2 — k 2 )/2 = a' 2 /2. Moreover, 
for | a' | < V22n/L, one finds kq 2 ~ q 3 . Then, the linearized equation ( 1B5I) becomes 



J2 



y y 



f. = Tf^m ■ (B6) 

Berger (1985) goes further on to derive a linearized expression for the total relative 
magnetic helicity in the volume V above S, namely 



Although not explicitly mentioned in Berger's (1985) analysis, equation ( 1B7I) appears to 
occur by assuming that q 2 — kq ~ a' 2 , instead of q 2 — kq ~ a' 2 /2, that we have assumed 
in equation (1B6|) . In our formulation, therefore, Berger's (1985) equation (1B7I) is a factor of 
two too high. If no linearization is performed, then our equation ( ]B5j) is in agreement with 
Berger's (1985) analysis. 

In an observed vector magnetogram, the best a-value is inferred in physical units of 
inverse length. The scaled value a' of the force-free parameter relates to a via the equation 
a' = ad where d is the elementary length in the magnetogram, expressed in physical units. As 
in any discrete parameter distribution with a well-defined (preferably fixed) length element, 
the length d in the magnetogram can be naturally represented by the linear size of the 
magnetogram's pixel. From this understanding and using the definition q 2 = u 2 + v 2 , the 
linearized expression fi for / can be written as 

\h 2 

EV~~* \ u u,v\ 

fi = Tid a where T\ = ,, 2 — . (B8) 

It is important to emphasize that the linearization fi as shown in equation ( 1B8I) provides 
a lower limit of /, and hence a lower limit of the free magnetic energy E c and the relative 
magnetic helicity H m in the LFF approximation (equations (35) - (37)). Green et al. (2002) 
and Demoulin et al. (2002) reached the same conclusion when deriving the linearized helicity 
expression of equation (J2J). Since the LFF magnetic energy is the minimum energy for a given 
relative helicity (a consequence of the Woltjer- Taylor theorem), the linearized energy/helicity 
expressions are underestimations of the actual energy/helicity values. The underestimation 
of E c and H m is negligible for small values of \a'\ and increases as \a'\ — > (2n/L). This, 
however, does not invalidate the linearized energy/helicity expressions for larger \a'\. As 
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we see in Figure [Tb, the underestimation factor is reasonable even for large \a'\, at least in 
view of other sources of uncertainties that are expected for observed magnetic configurations, 
and especially the use of the LFF approximation itself. Demoulin (2006), also provides a 
practical explanation of the underestimation effect based on well-known properties of the 
LFF magnetic fields. 
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Fig. 1. — Comparison between the volume- and the surface-integral expressions of the relative 
magnetic helicity H m in three LFF dipole models with analytical boundary conditions for the 
vertical magnetic field component. Boundary conditions are taken by the submerged poles 
(SP) model solution, the Sakurai-Uchida (SU) solenoidal model solution, and the Gold- 
Hoyle (GH) model solution (see §4.3 for details). Different model results are represented 
by different colors, (a) Dashed curves and triangles refer to results of the volume-integral 
helicity, equation (j28|) . while the exact surface-integral results, equation (40a), are shown 
by solid curves and rectangles. The linearized surface- integral results, equation (40b), are 
shown by dotted curves, (b) Dashed curves refer to the ratios between the volume-integral 
helicity and the linearized surface-integral helicity, while solid curves show the ratios between 
the exact surface-integral helicity and the linearized surface-integral helicity for each of the 
models. In all cases, the footpoint separation length L sep of the dipoles has been kept fixed. 
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Fig. 2. — Parts of the magnetic configuration of the two studied solar ARs. Shown are the 
heliographic magnetic field components of the ARs on the heliographic plane, (a) Disam- 
biguated photospheric vector magnetogram of NOAA AR 8844 as obtained by the IVM on 
2000 January 25 at 19:02 UT. A vector length equal to the tick mark separation corresponds 
to a horizontal magnetic field of 2300 G. (b) Disambiguated photospheric vector magne- 
togram of NOAA AR 9165 as obtained by the IVM on 2000 September 15 at 17:48 UT. A 
vector length equal to the tick mark separation corresponds to a horizontal magnetic field 
of 1760 G. Tic mark separation in both images is 10". North is up; west is to the right. 
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Fig. 3. — Timeseries of the magnetic flux and flux imbalance in the two studied ARs. 
Negative- (positive-) polarity fluxes are indicated by dashed (solid) curves, with readings on 
the left ordinate, while the red line corresponds to the relative magnetic flux imbalance in 
the ARs, with readings on the right ordinate, (a) NOAA AR 8844 (b) NOAA AR 9165. 
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Fig. 4. — Example calculation of a unique a-value in NOAA AR 8844 for a vector magne- 
togram obtained at 18:18 UT on 2000 January 25. (a) Various a- values obtained for various 
significance thresholds a (see text). The straight solid line indicates the flux-weighted aver- 
age of these a-values and the shaded area indicates the uncertainty in the calculation of this 
average, (b) Estimates of the total unsigned magnetic flux in the AR for various significance 
thresholds a. A threshold of la corresponds to a vertical magnetic field of 100 G and a 
horizontal magnetic field of 200 G. 
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Fig. 5. — Timeseries of the constant force-free parameter a in the two studied ARs. The 
dashed line and the surrounding shaded area correspond to the estimated average a-value 
and its uncertainties, respectively, (a) NOAA AR 8844 (b) NOAA AR 9165. 
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Fig. 6. — Magnetic energy and helicity budgets in NOAA AR 8844. Red (blue) curves 
correspond to the linearized (exact) surface-integral expressions. The error bars have been 
calculated from the linearized expressions, (a) Timeseries of the total relative magnetic 
helicity H m . The dashed line and the surrounding shaded area correspond to the linearized 
average value and its uncertainties, respectively, (b) Timeseries of the magnetic energy 
budgets in the AR. The potential magnetic energy is shown by the green curve. The total 
(free) energy and its uncertainties are shown by the solid (dashed) curves. 
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Fig. 7.— Same as Figure El but for NOAA AR 9165. 
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Fig. 8. — Comparison between the potential and the total magnetic energy estimates pro- 
vided by our analysis (equations ( TT4l) and (38a), (38b)) and the Virial theorem (equation (1221) 
for the two studied ARs. Shown with blue (red) curves are the potential (total) magnetic 
energy estimates. The linearized (exact) estimates for the total energy are represented by 
solid curves and rectangles (dotted curves). The Virial-theorem estimates are represented 
by dashed curves and triangles. The error bars correspond to the linearized expression for 
the total energy. For clarity in comparing the different energy values, we have applied a 
different scaling for the total energy (with readings on the left ordinate) than the scaling for 
the potential energy (with readings on the right ordinate), (a) NOAA AR 8844. (b) NOAA 
AR 9165. 
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NOAA AR 


$ (xlO 21 Mx) 


a (Mm- 1 ) 


H m (xlO* 
Exact 


12 Mx 2 ) 
Linearized 


8844 
9165 


5.1 ±0.2 
17.1 ±0.8 


0.023 ±0.006 
-0.024 ±0.006 


1.73 ±0.8 
-13.5 ±7.8 


1.48 ±0.4 
-12.8 ±3.7 


Ratio 


3.4 ±0.06 


1 


7.8 ±0.7 


8.6 ±0.5 



Table 1: Synopsis of the average magnetic flux, a- value, and relative magnetic helicity bud- 
gets for NOAA ARs 8844 and 9165. The third row refers to the ratio I-P9165/-P8844I between a 
given parameter -P9165 of NOAA AR 9165 and the respective parameter -P8844 of NOAA AR 
8844. 



NOAA AR 


E p (xlO 32 erg) 


E c (xlO 
Exact 


32 

Linearized 


E (xlO 
Exact 


32 

Linearized 


8844 
9165 


2.97 ±0.1 
9.52 ±0.5 


0.19 ±0.13 
1.44 ± 1.1 


0.15 ±0.1 
1.33 ±0.6 


3.16 ±0.2 
10.96 ± 1.2 


3.13 ±0.2 
10.85 ±0.9 


Ratio 


3.2 ±0.06 


7.6 ± 1 


8.9 ±0.8 


3.5 ±0.2 


3.5 ±0.1 



Table 2: Synopsis of the average potential, free, and total magnetic energy budgets, respec- 
tively, for NOAA ARs 8844 and 9165. The third row refers to the ratio (-P9165/-P8844) between 
a given parameter -P9165 of NOAA AR 9165 and the respective parameter -Ps844 of NOAA 
AR 8844. 



